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>• ' Abstract 

' This paper complements a new lattice formulation of SU(2) Yang-Mills the- 

l/^ . ory written in terms of new variables in a compact form proposed in the previous 

(N ■ paper. The new variables used in the formulation were once called the Cho- 

' Faddeev-Niemi or Cho-Faddeev-Niemi-Shabanov decomposition. Our formu- 

lation enables us to explain the infrared "Abelian" dominance, in addition 
! to magnetic monopole dominance shown in the previous paper, in the gauge 

invariant way without relying on the specific gauge fixing called the maximal 
Abelian gauge used in the conventional investigations. In this paper, especially, 
we demonstrate by numerical simulations that gluon degrees of freedom other 
I than the "Abelian" part acquire the mass to be decoupled in the low-energy 

region leading to the infrared Abelian dominance. 
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1 Introduction 



The dual superconductivity in Yang-Mills vacuum [T] is believed to be a promising 
mechanism j2] for quark confinement. For this picture to be valid, it is supposed that 
one can extract the relevant Abelian component responsible for duality from non- 
Abelian gauge theory, i.e., Yang- Mills theory, since the duality is understood as the 
electric-magnetic dual in the Abelian gauge theory represented by the Maxwell theory. 
Therefore, it is important to give a prescription to extract a variable playing the role 
of such an Abelian part from the original non-Abelian gauge theory in some sense. 
However, such a variable is not necessarily the Abelian gauge field in its appearance. 
To emphasize this situation, we use double quotation marks as "Abelian" . 

In the previous paper [3], we have proposed a formulation of Yang-Mills theory on 
a lattice in a compact form (say, compact lattice formulation). The compact lattice 
formulation was constructed so that it reduces in the (naive) continuum limit to the 
formulation of Yang-Mills theory written in terms of new variables obtained through 
non-linear change of variables (NLCV) from the original gauge field [1]. The NLCV 
generates the new variables which have the same form as the Cho-Faddeev-Niemi 
(CFN) or Cho-Faddeev-Niemi-Shabanov (CFNS) decomposition [SIEIE]- 

Prior to the compact lattice formulation, we have already given another lattice 
formulation in a non-compact form (referred to as the non-compact lattice formula- 
tion hereafter) [8]. The non-compact lattice formulation enabled one to define the 
magnetic monopole in the gauge invariant way in Yang-Mills theory on a lattice 
without introducing fundamental scalar fields. This is a remarkable result, since the 
conventional approach of defining the magnetic monopole in Yang-Mills theory with- 
out fundamental scalar fields heavily relies on a specific choice of gauge fixing, the 
so-called the maximal Abelian gauge (MAG) which breaks the color symmetry 
explicitly in addition to the local gauge symmetry. The gauge-invariant magnetic 
monopole in pure Yang-Mills theory has been constructed by introducing the unit 
vector field based on NLCV which plays the role of recovering color symmetry broken 
by a specific choice of color direction according to a Cartan decomposition in MAG. 
However, the magnetic charge resulting from the magnetic monopole defined in this 
way is not integer-valued in the non-compact formulation. 

This drawback was remedied by the subsequent compact formulation [3] which 
guarantees that the magnetic charge is integer-valued and obeys the Dirac quantisa- 
tion condition. Moreover, the infrared "Abelian" dominance and magnetic monopole 
dominance in the string tension were demonstrated by numerical simulations in the 
compact formulation, although such phenomena were found for the first time in the 
MAG [ini[IIl[l2l[13] . These results strongly support the dual superconductor picture 
of QCD vacuum as a promising mechanism of quark confinement. However, what is 
the mechanism for the infrared "Abelian" dominance or magnetic monopole domi- 
nance is an unanswered question in this investigation. From a theoretical point of 
view [13], on the other hand, it was clarified which variables should be identified 
with the "Abelian" part which is responsible for quark confinement within the 
continuum formulation [3] so that the variable gives the dominant contribution to 
the string tension to be calculated from the Wilson loop average (infrared "Abelian" 
dominance). This means that the remaining variable = — decouples in the 
low-energy or long-distance region to become irrelevant for the string tension, once 
such an identification of the "Abelian" part is achieved. 
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Our continuum formulation of Yang-Mills theory allows one to introduce the mass 
term |M|^X^ for the remaining field without breaking the local gauge invariance. 
Therefore, the dynamical mass generation for is not prohibited in this formulation. 
If such gluon mass for X^ is generated, it could be a mechanism for infrared "Abelian" 
dominance as pointed out in [11]. The dynamical generation of the gluon mass for X^ 
yields the decoupling of these degrees of freedom in the low-energy region leaving the 
"Abelian" part as the low-energy modes relevant to quark confinement. To confirm 
this scenario for dynamical Abelian projection is a main motivation of this paper. For 
MAG, mass generation was so far reported for the off-diagonal gluon component [T5] . 
In this paper, we have measured the mass for the gluon X^ directly by numerical 
simulations based on our lattice formulation [3]. In this paper, moreover, we discuss 
in detail how to define the lattice variable corresponding especially to the remaining 
part X^ in the compact lattice formulation to obtain the continuum counterparts to 
the errors of lattice spacing e. Preliminary results have already been reported in [16] . 
Finally, it is worth mentioning that the gauge invariance of the Abelian confinement 
mechanism was also discussed recently in many publications by other group, e.g., [17]. 



2 Compact formulation 

In order to consider a compact lattice version for the new formulation of Yang-Mills 
theory, we recall the continuum formulation presented in [1]. In the continuum for- 
mulation [3,11], we have introduced a color vector field n{x) = (nA(x)) {A = 1,2,3) 
of a unit length, i.e., n{x) ■ n{x) := nyi(x)n^(x) = 1. In what follows, we use the 
arrow to denote the vector and use the boldface letter to express the Lie-algebra 
SM(2)-valued field, e.g., n(x) := nA{x)TA with generators of sm(2), T4 = \aA where 
a A {A = 1, 2, 3) are Pauli matrices. The sM(2)-valued gluon field (gauge potential) 
A^(x) is decomposed into two parts: 

A,{x)=Y,{x) + X^{x), (1) 

in such a way that the color vector field n(x) is covariantly constant in the background 
field V^(x): 

= %[V]n{x) := 9^n(x) - ig[V^{x), n(x)], (2) 

and that the remaining field X^(x) is perpendicular to n(a;): 

= n{x) ■ Xf,{x) = 2tr(n(x)X^(x)). (3) 

Here we have introduced the gauge coupling g and we have adopted the normalization 
for generators: tr(T4Te) = ^Sab- Note that n'^(x), V^(x), X^(x) and A^(a:) are 
real-valued fields and their Lie-algebra forms are Hermitian due to Hermiticity of the 
generators T^. 

By solving the defining equation (I2|), the V^(x) field is obtained in the form: 

V,{x) = Vil(x) + V;^(x) = c,ixMx)-tg-'[d,nix),iiix)], (4) 

where the second term V^(a:) := — i(7~^[9^n(x), n(x)] = g~^{d^n{x) x n{x))ATA is 
perpendicular to n(x), i.e., n{x) ■ V^{x) = 2tr(n(x)V-[(x)) = 0. Here it should 
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be remarked that the parallel part V|[(x) = c^(a;)n(a;), c^{x) = tr(n(x)A^(x)) pro- 
portional to n(x) can not be determined uniquely from the defining equation 
Imposing the perpendicular condition determines V|l (x) and the remaining part 
X^(x) as 

X^(x) = -ig-'Hx), %[A]nix)]. (5) 

It is easy to check that the sum of V^(a;) and X^(a;) specified respectively by and 
dS]) agrees with the original field A^(a;) according to ([T]). 

On a lattice, on the other hand, we introduce the site variable n^. constructed 
according to [8], in addition to the original link variable L^^.^t- Note that we define a 
color vector field n^; := u^cta on the lattice corresponding to the continuum notation 
n{x) := nA{x)TA. In this paper, we define the link variable Ux,fj, as the exponential 
of the line integral of a gauge potential A^(x) along a link from x to x + fi: 

U,^^ = ^ exp (^-tg ^"^"^ dx^A^{x)J , (6) 

where e denotes the lattice spacing and ^ denotes the path ordering. In the ex- 
plicit estimation of the naive continuum limit, we adopt in this paper the mid-point 
definition for the link variable: 

Ux,f, = exp{-igeA^{x')), (7) 

using the midpoint x' := {x + efi/2, fi) of the link [x, x + e/i) running from x to x -|- e/i. 
This prescription is adopted to suppress as much as possible lattice artifacts coming 
from a finite (nonzero) lattice spacing, in contrast to our previous paper [3] where we 
have adopted the very naive definition: Ux,^ = exp{—igeA^{x)) . 

The link variable Ux^n and the site variable rix transform under the gauge trans- 
formation II [1] as 

Ux,^, ^xUx,^,^l+^, = U'^,^, ^ ^x^^x^l = K- (§) 

Note that is Hermitian, nj. = n^, and f/x,^ is unitary, t/J^^ = U^J^. It should 
be remarked that this transformation property follows from the most general form 
([6]) for the link variable U^^fj,, irrespective of the prescription for the discrete lattice 
approximation. 

The lattice variables V^^^ and X^^fj, corresponding to V^(x) and X^(x) should be 
expressed in terms of the site variable n^, and the original link variable U^^^, just as 
the continuum variables V^(x) and X^(x) are expressed in terms of n(x) and A^(x), 
However, the definition of lattice variables Vx^/j, and Xx^fj, is not unique. They must 
be defined in a consistent way with the defining equation on a lattice respecting the 
transformation property. We achieve this by solving a lattice version [3] of ([2]) and 



nxK,M - Vx,f,nx+^, = 0, (9) 

tr(n,X,,^) = 0. (10) 

The defining equation must be invariant under the gauge transformation, namely, 
they are form-invariant: n^^'^^ — Vx,^J.^'x+|J, — 0' ^^('^x^x,^) — 0- identify 
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the lattice variable \4^^ with a link variable which transforms in the same way as the 
original link variable U^,^'- 

This requirement guarantees that the defining equation ([9]) is gauge invariant. On 
the other hand, we define the lattice variable X^,^ so that it transforms in just the 
same way as the site variable n^;: 

Xx,fj. VLr^Xx^^VL^^ = X'^^^, (12) 

to realize the adjoint color rotation at the site suggested from the transformation 
property of the continuum variable. By this choice, indeed, the orthogonality condi- 
tion (fTOj) is kept gauge invariant. 

Explicit construction of the new lattice variables are as follows. We define V^^^ as 
a link variable which is a group element of G = SU{2) related to the sM(2)-valued 
background field V^(x) through 

K,^ = ^ exp (^-ig ^^''^ dx^V^ix)^ . (13) 

In the mid-point definition the link variable Vx^fj, reads 

Vx,f, = exp{-iegY^{x')), (14) 

where V^(x') is to be identified with the continuum variable V^(x) defined by (jlj) in 
the continuum limit. Hence K,.^^ must be unitary ^ = V~^. The same remark as the 
the link variable f/^,^ for the naive continuum limit holds also for the link variable Vx,^- 
In the previous paper [3] , the lattice version ([9]) of the defining equation ([2]) has been 
solved and the resulting link variable K;,^ is of the form (up to the normalization) [18] : 

K,/. = Vx,t,[U, n] = Ux,t, + n^t/^,/,n^+^, (15) 
and the unitary link variable Vx^^\U, n] has been obtained after the normalization: 



= Vx,^[U, n] := VxJ^iT[Vl,Vx,^]l2. (16) 



Indeed, the naive continuum limit e of the link variable (1161) reduces to the 
continuum expression (jlj). 

A naive choice for the lattice variable X^^^ is given by U^^/jV^^^ or V^-^^^Ux-^,^- 
These are suggested from the relation X^(x) = A^j(x) — V^(x) = —Y^{x) + A^(x). 
In fact, they satisfy the desired transformation property (fT2l) . Note that V^JJx,^i or 
Ux~fi,^iVx-ij,^ij_ is excluded, since it obeys the adjoint rotation at x + /i or x — /i, not at 
X. Then we can construct a lattice variable Xx^/j, as the linear combination: 

Xx,^ = aUx,^yl, + PVl^^^Ux-^.,^., (17) 



to satisfy the desired transformation property (1121) . Now we can see that it is rea- 
sonable to adopt ( ITOl) as a lattice version of the orthogonality equation ([3]). In fact, 
tr(na;X^(x)) = implies (See Appendix for the derivation of (ITSl) ) 

tr(n,X,,^) =(« + /5)tr(n,{l - ige^^{x)}) + (« - p)0{e^) 

= - ig{a + /5)etr(n^.X^(x)) + (a - p)0{e^) = + (a - I3)0{e^). (18) 
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In this way, the lattice variables X^^fj, is expressed in terms of the site variable and 
the original link variable U^^^ as 

= Xr^^^lU, n] = {a + P)l + aUr,^^nr,+f,Ul^n^ + f3nJjl_^^^\i^_^Ux-^^^. (19) 

In particular, a good choice is obtained for the symmetric case, i.e., a = /3, since 
this choice enables us to define the lattice variable X^,^ so as to reproduce the naive 
continuum limit of the orthogonality equation up to 0{e^). 

Finally, we obtain the unitary lattice variable X^,^[f/, n] after the normalization: 

X,,^[f/, n] :=X,,^/v/tr[Xl,^X,,^]/2. (20) 

In the numerical simulations, we have adopted the choice a = (3 hy the reason men- 
tioned above. Moreover, there are some arbitrariness for extracting the Lie-algebra 
valued variable from the compact lattice variable Xa;^^[?7, n]. This issue will be 
examined by comparing the results of numerical simulations. 

3 Abelian dominance and gluon mass generation 

3.1 Identifying the Abelian part and mass term for the re- 
maining part 

The field can be regarded as the "Abelian" part in the reformulated Yang-Mills 
theory by the following reasons. 

i) The "Abelian" part corresponds to the diagonal part of the gauge potential 

in the context of the conventional MAG which is reproduced when the color 
vector is aligned in the same direction over the whole spacetime, for example, 

n(a:) ^ no := (0,0,1). (21) 

ii) The Wilson loop average W{C) in Yang-Mills theory written in terms of A^ is 

rewritten into the reduced Wilson loop average W{C) which is entirely rewritten 
in terms of in the reformulated Yang-Mills theory, as demonstrated in [H] . 

iii) The mass term for can be introduced without breaking gauge invariance in 
this reformulation [3]. In fact, it has been shown to one- loop order [13] that 
such an effective mass term is generated due to the gauge-invariant dimension 
two condensate <^Xj^X^^ thanks to the gauge invariant self-interaction term 
i^2^gABC^B^C'^2 a,mong X^ gluons, in sharp contrast to the ordinary self- 
interaction term ^g"^ [e^^'~^ A.^ A!^)"^ which is not gauge- invariant. 

Therefore, in the energy region lower than the mass Mx of the field X^, the remaining 
components X^ should decouple or negligible and the field could be dominant. 
This leads to the infrared Abelian dominance (in the string tension) in our reformu- 
lation. 
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Keeping these facts in mind, we proceed to obtain a fitting function of the two- 
point correlation function. Suppose that the Yang-Mills theory has the effective mass 
term: 

^ r2 vA^A (^22) 



An additional quadratic term in of the following type could be generated from 
gauge fixing conditions in the differential form [Ij. 



2/3 



A\2 



(23) 



This can be understood as follows. Recall that we impose an constraint called the 
new Maximal Abelian gauge (nMAG) to obtain the reformulated Yang-Mills theory 
with the original gauge symmetry SU(2) even after introducing the color field n(x) 
which apparently increases gauge degrees of freedom [1] . Then we introduce a gauge- 
fixing parameter a for nMAG of the form: — ^^((^^[VJX^)'^)^. This term does not 
fix the SU(2) gauge invariance. Therefore, we adopt the Landau gauge for the overall 
gauge fixing of whose differential form is d^A.^ = 0. This gives an additional 
quadratic term: —^{d'^X.^)'^ coming from the GF term: —^{d^A^)'^. Therefore, 
combining two terms yields an additional term quadratic in X^: —^{d^X^Y with 
f3~^ = a"^ + a'~^. Thus we assume the effective propagator for X gluon of the form: 



jjXX 



{k) 



Ml 



'IIP 



kfiky 



k^ - pM[ 



(24) 



In particular, the limit /5 — + oo reproduces the Proca case: 



k^ - Ml 



e k^k 



Ml 



(25) 



This form was adopted in the study of off-diagonal gluon mass generation in MAG [T5] 
where the mass term iMl^A^A^ was introduced by hand without preserving the gauge 
invariance. Note that both nMAG and Landau gauge conditions are exactly satisfied 
only at a = and a' = 0. This is realized at /3 = limit: 



r.XX 



(k) 



k^ - Ml 



kfj,ki^ 



k^ 



rjXX 



{k) 



A;2 - Ml 



D 



XX 

flflOO 



{k) 



Ml 



■ (26) 



Therefore, the /3 = limit differs from the previous Proca case used in MAG. However, 
it will turn out below that the constant shift of the propagator gives the same decay 
rate and hence the same mass Mx of X^ gluon. 



3.2 Numerical simulations 

We have generated configurations of link variables {U^,^} based on the standard heat 
bath method for the standard Wilson action. The numerical simulation are performed 
at /3 = 2.3, 2.4 on 24^ lattice, aX (3 = 2.3, 2.4, 2.5 on 32^ lattice, at /5 = 2.4, 2.5, 2.6 
on 36^ lattice, and at /3 = 2.4, 2.5, 2.6 on 48"* lattice by thermalizing 15000 sweeps. 



6 



r[fnn] 

0.2 0.4 0.6 




r[fm] 

0.6 0.8 



24b2.3J 
■32b2.4 
36b2.5 
43b2.6 
■24b2.3 
32b2.4 
36b2.5 
■48b2.6 
.[24b2.3 
) 32b2.4 
: 36b2.5 
) 48b2.6 
'24b2.3 
■32b2.4 
36b2.5 
48b2.6 



0.1 

■go.01 

0.001 
1e-04 



■i 




0.5 1 1.5 



Figure 1: Logarithmic plots of scalar-type two-point correlation functions Doo'i'f') '■= 
{0{x)0'{y)) as a function of the Euclidean distance r := \J{x — yY for O and O' . 
(Left panel) 0{x)0'{y) = V;t(x)V^(y), A;^(x)A;^(y), -V;^(x)X;^(2/), X;^(a;)X^(|/), 
(Right panel) C(x)0'(?/) = n'^(x)n^(?/), c^{x)c^{y), X^(x)X^(y), from above to below- 
using data on the 24^ lattice (/3 = 2.3,2.4), 32^ lattice (/5 = 2.3,2.4), 36^ lattice 
{P = 2.4, 2.5), and 48^ lattice (/? = 2.4, 2.5, 2.6). Here plots are given in the physical 
unit [fm] or in unit of square root of the string tension ^ciphys- 



Here 200 configurations are stored every 300 sweeps. Other settings of numerical 
simulations are the same as those in the previous paper [3]. 

We are now ready to study characteristic features of the reformulated Yang-Mills 
theory written in terms of new variables n'^(x), c^{x), X^(x) defined through NLCV 
of the original field variable A^{x): infrared Abelian dominance, magnetic monopole 
dominance and non-vanishing gluon mass. Among them, the magnetic monopole 
dominance in the string tension has already been confirmed in the previous paper [3] 
using the gauge-invariant magnetic monopole which is guaranteed to have integer- 
valued magnetic charge subject to the Dirac quantization condition according to our 
construction of magnetic current based on NLCV. An advantage of our formulation 
is that we can confirm such characteristic features for any choice of gauge fixing, not 
restricted to MAG, since our formulation allows us to take arbitrary type of gauge 
fixing for the original variable A^(x). 

To study the infrared Abelian dominance and the gluon mass generation in the re- 
formulated Yang-Mills theory, we first define the two-point correlation functions (full 
propagators) for the independent variables in the new formulation on a lattice, i.e., 
Cx,^i and X;^^, in addition to the original variable A;^^. For simplicity, we exam- 
ine just the contracted scalar-type propagator simplified by avoiding the complicated 
tensor structure: 

Dnn{x -y) = (n^ n;^) , Dcc{x - y) = {Cx',^ Cy,^^) , 

Dxx{x -y)= (Xj^ X^^) , Dx'x'ix - y) = (xf, ^ X^, ^) , (27) 

and 

^AA(x-y) = (Af,^ A^,^). (28) 

^ The proposed NLCV enables one to extract the "Abehan part" V;^^ irrespective of the choice 
of the gauge fixing preserving the color symmetry. The Yang-Mills theory in the conventional MAG 
is reproduced as a very special limit (PT|) of our reformulated Yang-Mills theory based on NLCV. 
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Figure 2: 



Logarithmic plots of the rescaled correlation function r^^'^Dooi'^) a-s a 



function of r for O = V^, Aj^, c^, (and X'^) from above to below, using the same 



colors and symbols as those in Fig. [H Here two sets of data for the correlation function 
Dxx{x — y) are plotted according to the two definitions ([30 
on a lattice. 



and (EB of the X;^ field 



Here the Lie-algebra valued gauge potential A^,/^^ or V^;'^^ is defined from the respective 
link variable by 



{i/2ge) 



a. 



ul 



X,fJ, 



{t/2ge) 



V, 



X,fl 



(29) 



For the variable X^^^, on the other hand, we examined two options: one is extracted 
from decomposing the gauge potential (group- valued) : 



{i/2ge) 



and the other is from the definition of the decomposition (Lie-algebra- valued): 



^x',fi 



^X' ,fl 



^x',^^. 



The field Cx\fj, is defined by 



^x' ,IJ, 



tr(n^K',M) = tr(K-,Mn 



^'-'■x+iij 



(30) 



(31) 



(32) 



The numerical results are presented in Fig. [TJ As is quickly observed from the 
left panel of Fig. [1], Dvv{x — y) and Daa{x — y) exhibit quite similar behaviors in 

In order 



\x - y\ 



the measured range of the Euclidean distance r 
to determine the physical scale, we have used the relationship between the (inverse) 
gauge coupling /3 and lattice spacing e given in [11] El which is summarized in Table [H 
From the right panel of Fig. [T], Dyvi^ ~ v) {Dnn{x — y) or Dcc{x — y)) is dominant 
compared to Dxx{x—y) which decreases more rapidly than other correlation functions 
in r. This implies the infrared "Abelian" dominance, provided that the components 
V^(x) composed of n;^ and Cx^^ are identified with the "Abelian" part of A^{x). 



^ Wc use the relationship between the physical units, IGeV 
5.06773fm"^ This comes from he = 0.197327GeV • fm. 



0.197327fm or IGeV 
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lattice spacing e 


lattice size L [fm 




p 


[1/ V^p?*?/*] 


[fm] 


24^ 


32^ 


36^ 


48* 


2.3 


0.35887 


0.1609 


3.8626 


5.1501 


5.7939 


7.7252 


2.4 


0.26784 


0.1201 


2.8828 


3.8438 


4.3242 


5.7657 


2.5 


0.18551 


0.08320 


1.9967 


2.6622 


2.9950 


3.9934 


2.6 


0.13455 


0.06034 


1.4482 


1.9309 


2.1723 


2.8964 



Table 1: The lattice spacing e and the lattice size L of the lattice volume at 
various value of (3 in the physical unit [fm] and the unit given by ^JOphys- 

As is seen from the left panel of Fig. [1], a non-trivial mixed correlation function 
V^(x')X^(y')^ < exists, since V^(a;) includes a perpendicular component to w^{x). 
Fig. [T] demonstrates nice independence of our results against variations of the ul- 
traviolet cutoff (the lattice spacing e). The propagators calculated at the lattices with 
different e follow the same curve if plotted in the physical units. These accurate plots 
provide an additional support that the results presented in this paper are definitely 
not lattice artifacts. 

Note that we must impose the gauge fixing condition for the original variable Aj;/^^ 
to obtain the correlation function. In our simulations, we have chosen the lattice 
Landau gauge (LLG) for the original field A^(x) for this purpose. Thus we have 
confirmed the infrared "Abelian" dominance with color symmetry being kept, since 
the Landau gauge keeps the color symmetry. This is one of our main results. The 
infrared Abelian dominance was so far obtained only for the MAG which breaks the 
color symmetry explicitly. As already mentioned, moreover, we can choose any other 
gauge and we can study using this formulation if the infrared "Abelian" dominance 
can be observed in any other gauge. We hope we can report the results in the other 
gauge in future investigations. 

Next, we determine the gluon mass generated in the non-perturbative way by 
examining the correlation functions in more detail. The gauge boson propagator 
D-^^[x — y) is related to the Fourier transform of the massive propagator D-^^[k\. 

D^J'ir) = {x^^ix)X^iy)) = J ^^.e^'^-'^D^J'ik). (33) 

Then the scalar-type propagator D^-^{r) := D-^^i^x) as a function of r should behave 
for large Mxr as (See [15] for details of the integral calculation.) 



D..(r) = (x-(x)X-(,)) = 7^e-(^-)-^ . (34) 



(27r)4^ P + ~ 2(27r)3/2 r^/s 

Therefore, the scaled propagator r'^l'^Dxxij'^ should be proportional to exp(— M^r) 
for MxT ^ 1 with M being the damping rate of r'^l'^Dxx{j'\ In other words, the 
mass Mx of the gauge field can be estimated from the slope in the logarithmic 
plot of the scaled propagator r'^l'^Dxxij'^ as a function of r. E| 



^ Here we have assumed that the anomalous dimension is sufficiently small so that the exponent 
of the power of r is the same as the tree value. 
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Figure 3: Gluon "mass" and decay rates (in units of GeV and y/a^h^) as the 
function of the inverse lattice volume 1/V in the physical unit. (Left panel) for 

from above to below extracted according to the fitting: 
exp(— M^r), (Right panel) for n'^(x) extracted according to the 



-3/2 



{0{x)0{y)) 
fitting: (n^(x)n'^{y)) ~ exp(— M„r). 



Fig. [2] shows the logarithmic plots of the scaled scalar-type propagator for A^',^, 
Cx'^fi and lLx,ti as a function of the distance r measured in the physical unit [fm] 
and in unit of square root of the string tension ^cXphys = 440 MeV. According to 
Fig. [21 we find just small difference between two types of Dxx{x — y) defined by 
(!30|) or fl3T!) over several choices of lattice spacing (i.e., several values of P, P = 
2.3,2.4,2.5,2.6). Therefore, we can use either definition of the lattice variable X^,/^^ 
to obtain Dx'x'i^ — y) in the consistent manner. 

In Fig. [H], the measured values for the gluon mass are plotted as the function of 
the inverse lattice volume 1/V in the physical unit, to study the finite-size effect on 
the mass. The finite lattice-size effect seems to be small for the gluon mass Mx- Here 
the error bars originate from the fitting procedure for obtaining the slope, but no 
systematic errors such as finite-volume are included. In this way, we have estimated 
the mass for the X gluon: 

Mx ^ 2.98^0^5;^ ~ 1.31GeV, 

Mx' ^ 2.69^/ap^ ~ l.lQGeV. (35) 

Even after the whole gauge fixing, our formulation preserves color symmetry in 
sharp contrast to the conventional MA gauge. In view of the fact that our reformu- 
lated Yang-Mills theory reproduces the Yang-Mills theory in MA gauge as a special 
limit, the remaining part X^(x) could correspond to the off-diagonal part in this 
limit. From this point of view, our result is consistent with the result obtained for 
the off-diagonal gluon mass in MAG [T5] . 

Moreover, we have simultaneously estimated the decay rate for the new fields 
n'^(x), c^(x), Vj^(x) and the original gauge field A^{x) by imposing the LLG as 
the overall gauge fixing. For O = X^,c^,A^,V^, the decay rate Mq is extracted 
according to the fitting: {0{x)0{y)) ~ r~^/^ exp(— M^r). Fig. [S] indicates not so 
small finite volume effect for data of 1/V > 0.02. Using the data of 1/V < 0.02, 
therefore, we have estimated the decay rate (or "mass" ) as 

M„ ~ 2.24 y^pi^ ~ 0.986GeV, 
Me ~ 1.94y5pJ^ ~ 0.856GeV, 

Ma ~ 1.35y/a~^ - 0.596GeV. (36) 



10 



The decay rate obtained from the correlation function of c^{x) field is slightly 
larger than that expected from the result in MAG. It should be remarked that the 
decay rate for the correlation function of n"^(x) field is extracted according to the 
fitting function (n"^ {x)n"^ (y)'^ ~ exp(— M„r) which is not yet justified from the the- 
oretical consideration. This might be an origin of the large value of M„. More 
simulations on the larger lattice are expected to eliminate finite volume effect for 
these values. However, we have no argument for guaranteeing the gauge invariance 
of these values or for identifying these values with their "masses". In fact, the field 
Cn{x) is not gauge invariant. These issues will be checked in further investigations 
based on our reformulation. 

Finally, we comment on the "Abelian"part V^(x), since our treatment of the 
"Abelian"part V^(x) is different from the conventional approach based on MAG. The 
above result yields the "mass" of the "Abehan"part Y^{x) : My ^ Ma ^ 0.59GeV . 
This value is nearly equal to that of the diagonal gluon mass obtained by imposing 
the Landau gauge in the conventional approach as reported in the second paper of [I5] 
where the Landau gauge was imposed on the Abelian diagonal part a^(x) in addition 
to the MAG for off-diagonal gluon field defined by the Cartan decomposition 

A^{x) = A^(x)T" + a^{x)T^ (a = 1,2). Therefore, the prescription of gauge fixing 
in [15] is different from ours. 

4 Conclusion and discussion 

In this paper we have developed a compact lattice formulation of SU(2) Yang-Mills 
theory proposed in the previous paper [3] as the lattice version of the NLCV which 
was once called the CFN or CFNS decomposition. This resolves all drawbacks of the 
previous non-compact lattice formulation of our own [8]E This compact formulation 
has enabled one to define the gauge-invariant magnetic monopole with the magnetic 
charge subject to Dirac quantisation condition and to extract the "Abelian" part 
V^(x) yielding the infrared "Abelian" dominance in the string tension for any choice 
of the gauge fixing for the original gauge field A^(x) in the original YM theory. 

In order to confirm the dynamical mass generation for the remaining part X^(x) as 
a mechanism for the infrared "Abelian" dominance, we have measured the two-point 
correlation function (the full propagator in real space) in our lattice formulation by 
imposing LLG for the original gauge field A^(x) as the whole gauge fixing. We have 
found the infrared "Abelian" dominance in the sense that the X^(x) propagator is 
suppressed in the long distance compared to n(x) and c^(x) (and V^(x) propagators 
as an immediate consequence of dynamically generated mass Mx = 1.2 ~ 1.3 GeV 
for X^(x) (which is larger than the decay rate of other gluon field propagators). 

Even after the whole gauge fixing, our formulation can preserve color symmetry 
by choosing the gauge-fixing condition which does not break color symmetry, e.g.. 
Landau gauge. This opens a path to examine color confinement in the same frame- 
work as quark confinement in the dual superconductivity picture. This feature is 
in sharp contrast to the conventional MA gauge breaking color symmetry, although 
our formulation reproduces the MA gauge as a special limit (12T]) . It is important to 

^ This is done up to specifying the integration measure for the respective new variable. It was 
not necessary to resolve this issue for obtaining the results reported in this paper. 
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demonstrate explicitly the gauge-fixing independence of our results obtained in this 
paper for establishing the gauge-invariant mechanism for quark confinement. 
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A The accuracy of the naive continuum hmit 



For the naive continuum limit e — 0, we show that the lattice variable defined by 
:= + P^x-tM,f,U:,:-f,,^ yields X^;,^ = X^(x) + 0{e^) for a = P, while 



X 



X^(x) -|- 0{e) for a ^ f3. The repeated use of the Baker-Campbell-Hausdorff 



formula yields 

= exp{-igeAx+|,/2,|,)exp{ige¥x+^,/2,^,) or exp{ige¥x-^,/2,f,) exp{-igeAx-.f,/2,f,) 



exp <^ -igeXx±f,/2,f, ± 



exp < -igeXx^^ =F i—df^X^^^ ± 



(ge)" 



x^li/2,^1 ^ x±fi/2,ii 
2 



+ 0(6^) I 

x±^i/2,n,^x±^i/2,fi + O(e^) 



exp(-z(7eX,,^) exp j^z^^^X,,^ ± ^ [A,,^, + ^(e^) 
1 T ± [A.,^, + 0(e3) 



exY>{-igeKx,^) 



(37) 



where we have used X^^,^ := A^,^ - V^,/,, V^.±^/2,^ = V^^.,^ ± e/2d^Yx,^, Xx±^/2,^ 
Xx.,;t ± e/2dfjXx,^i, and A^.±^/2,^ = A^.,^ ± e/29^A^,^. Thus we have 



= exp{-igeXx,^ 



.ge' 



(gey 



{a + P) exp 



a + P)l + {a-P)^-t^d,Xx,,+ 2 



■ge^ 



[Ax.,^,V,j| + 0(e' 

(gey 



-tgeXx,, + {a- f3)/{a + f3) ^-z^d.Xx,, + ^ [A,,^, V,.j| + 0(e=^) 

(38) 



The statement follows from the fact that the choice a = P eliminates order terms. 
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